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Abstract 

A recently developed time dependent model for the excitation of a 
nucleon from a bound state to a continuum resonant state in the sys- 
tem n+core is applied to the study of the population of the low energy 
continuum of the unbound 10 Li system obtained from 11 Li fragmen- 
tation. Comparison of the model results to new data from the GSI 
laboratory suggests that the reaction mechanism is dominated by final 
state effects rather than by the sudden process, but for the population 
of the 1=0 virtual state, in which case the two mechanisms give almost 
identical results. There is also, for the first time, a clear evidence for 
the population of a d 5 / 2 resonance in 10 Li. 

1 Introduction 

In this paper we apply the projectile fragmentation model developped in Ref. 
[1] to the the elastic breakup (diffraction reaction) of the two-neutron halo 
nucleus n Li on a 12 C target, as recently measured at GSI and presented 
in Ref. [2]. The observable studied is the one- neutron-core relative energy 
spectrum. These type of data are necessary to enlighten the effect of the 
neutron final state interaction with the core of origin. 
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Light unbound nuclei have attracted much attention [5]- [JS] in connection 
with exotic halo nuclei. Besides, a precise understanding of unbound nuclei 
is essential to determine the position of the driplines in the nuclear mass 
chart. In two-neutron halo nuclei such as 6 He, n Li, 14 Be, the two neutron 
pair is bound, although weakly, due to the neutron-neutron pairing force, 
while each single extra neutron is unbound in the field of the core. In a 
three-body model these nuclei are described as a core plus two neutrons. 
The properties of the core plus one neutron system are essential to structure 
models which rely on the knowledge of single particle quantum numbers such 
as angular momentum, parity, energies, as well as the corresponding neutron- 
core effective potential. With those ingredients the spectroscopic strength 
function of neutron resonances in the field of the core can be obtained. 

We will use the projectile fragmentation formalism, an inelastic-like ex- 
citation to the neutron-core continuum, to the study of the effect of final 
state interaction of the neutron with the projectile core. The model is a 
theory which has already been shown to be relevant to the interpretation of 
neutron-core coincidence measurements in nuclear elastic breakup reactions 
with projectiles of 14 B and 14 Be [TJ [2]. In the case of two nucleon breakup we 
describe only the step in which a neutron is knocked out from the projectile 
by the neutron-target interaction to first order and then re-interacts in the 
final state with the core. The case in which a resonance is populated by a 
sudden process while the other neutron is stripped has been already discussed 
in Ref. [3] and it has been shown [lj that there is a simple link between the 
two methods of Ref. [1] and Ref. [3]. We assume that the neutron which is 
not detected has been stripped while the other suffers an elastic scattering 
on the target. The influence of the second nucleon is taken into account only 
by a modification of the neutron-core interaction in the final state. Section 2 
contains a brief reminder of the formalism from Ref. [TJ. Section 3 describes 
the results of our numerical calculations for 11 Li which, being already well 
understood [2]-[24j, is used here as a test case. Details on our assumptions 
for the potentials needed in the calculations are also presented. Finally our 
conclusions are contained in Sec. 4. 
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2 Formalism for inelastic excitation to the 
continuum. 

Following Refs. [H [501 EH [52], we will describe the inelastic-like excitations 
of one neutron, from a bound state ipi to a final state ipf in the continuum, 
by the time dependent perturbation amplitude : 



V2 is the interaction responsible for the neutron transition (cf. Eq. (2.15) 
of [Hn]). The potential V^r — R(£)) moves past on a constant velocity path 
with velocity v in the z-direction with an impact parameter b c in the x- 
direction in the plane y = 0. This assumption makes our semiclassical model 
valid at beam energies well above the Coulomb barrier. This is in fact the 
regime in which projectile fragmentation experiments are usually performed 
(cf. Sec. 3). The coordinate system used in the calculations is shown in 
Fig. [T] and it corresponds to the no-recoil approximation for the core. Let 
ipi(r,t) = (fii(r)e~% £it be the single particle initial state wave function. Its 
radial part 0«(r) is calculated in a potential Vws(r) (cf. Sec. 3.1) which 
is fixed in space. For initial and final states of different angular momentum 
our wave functions are trivially orthogonal due to the orthonormality of their 
angular parts. For transitions conserving the angular momentum of the single 
particle states the orthogonalization correction has been estimated in Ref. 
PQ and found negligible. Now change variables and put z — vt = z' or 




(1) 




target 



Figure 1: Coordinate system used in the calculations 
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t = (z — z')/v, define 

then choosing V^(r) to be a delta-function potential V 2 (r) = V26(x)5(y)S(z), 
with t> 2 = [MeV fm 3 ], the integrals over x and y can be calculated giving 

A ^ = r dz W« °> *)&( 6 « °> ^) el9Z - ( 3 ) 

iTlV J-oo ' 

The value of the strength t> 2 used in the calculation is discussed in Sec. 3. 

The following forms for the wave functions will be used. For the initial 
bound state 

<Mr) = -C^fW^^. # (4) 
Due of the strong core absorption implied in the following by Eq. (jSJ) , we use 
in this paper the asymptotic form of the initial state wave function, given in 
terms of the Hankel functions h\ ; . However the exact wave function, numer- 
ical solution of the bound state Schrodinger equation can be used without 
introducing further complexity in the calculations. For the final continuum 
state 
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Mr) = C f k-{h\-\kr) - S lf h%\kr))Y lf , m ,(e,<l>), (5) 

Cf = y2/L is the normalization constant for the final state. L is a large box 
radius used to normalize the continuum wave function (cf. Eq. (2.5) of Ref. 
|49j). Si f = e 2l5l i is the neutron-target S- matrix in the If partial wave. 

The probability to excite a final continuum state of energy e f is an average 
over the initial state and a sum over the final states. Thus introducing the 
density of final states and the density of final states, according to Ref. [19] 
p{ef)def = Lm / nh 2 kde f , the probability spectrum reads 



dP in _ 2 V 2 „ 2 m 1 |1_Q |2|r 12 (c\ 

de f - vrftV i n 2 k2i i + i m - mfl bm - m ^ 1 m ' m ^ ' w 

where now S = Se 2lu = e 2l ^ +u ^ is an off-the-energy shell S-matrix and v is 
the phase of the integral: 



Iirii,mf 



dze^h^^Y^M 0)k l -h\;\kr)Y lf , mf (9, 0). (7) 
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A detailed derivation of the above equations can be found in Sec. 2 of Ref. pQ. 
For simplicity the equations in this paper are given without spin variables 
in the initial and final states. The generalization including spin is given in 
Appendix B of Ref. [1] . 

The cross section differential in Ef is given as 



(see Eq. (2.3) of [53]) and C 2 S is the spectroscopic factor for the initial state. 

The core survival probability P c t(p c ) = \S ct \ 2 [53] in Eq. (jSJ) takes into 
account the peripheral nature of the reaction and naturally excludes the 
possibility of large overlaps between projectile and target. P ct is defined in 
terms of a S-matrix function of the core-target distance of closest approach 
b c . A simple parameterisation is P c t(b c ) = e^~ ln2exp ^ Ra ~ bc ^ a ^ [53], where the 
strong absorption radius R s ~ 1A(A^ 3 + A^ 3 ) fm is defined as the distance 
of closest approach for a trajectory that is 50% absorbed from the elastic 
channel and a=0.6 fm is a diffuseness parameter. 

3 Applications 

3.1 One neutron average potential 

We apply the fragmentation model to the study of the relative energy spec- 
trum n+ 9 Li obtained by the authors of Ref. [2] in the breakup reaction of 
n Li on 12 C at 264 A.MeV. The structure of n L i is already well known from 
a number of experiments and theoretical papers [2]- [21]: the two neutron 
separation energy is 0.3±0.27 MeV [31]; the wave function is combination of 
a 2s state with a spectroscopic factor 0.31, a pi/ 2 with spectroscopic factor 
0.45 [1] and there is also a small ds/2 component. The main d§/2 strength is 
in the continuum centered around 1.55 MeV [2]. The link between reaction 
theory and structure model is made by the neutron-core potential determin- 
ing the S-matrix in Eq. d3J. Then if the theory fits the position and shape 
of the continuum n-nucleus energy distribution, obtained for example by a 
coincidence measurement between the neutron and the core, the parameters 
of a model potential can be deduced. 

The calculations in the present paper are made with different potentials 
for the initial and final state. The initial wave function for the s-state is 
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calculated in a simple Woods-Saxon potential with R= tqA 1 ^ 3 and strength 
fitted to the separation energy 0.3 MeV[Sl] and whose parameters are given 
in Table [TJ 

V WS (r) = 1 + e(r . R)/a ~ I — J — il + e (r-R)/a ) 2 l - a ( 9 ) 



Table 1: Asymptotic normalization constants Ci(fm -1 / 2 ) for the initial state components 
wave functions of the bound neutron. Spectroscopic factors from Ref. 0]. Potential 
parameters are: Vo fitted to give the separation energy 0.3 MeV [53]. The other potential 
parameters are rp=1.27 fm, a=0.75 fm, V so =5.25 MeV. 





-0.3 MeV 


H 7 ji 




C 2 S 


1/2 


0.76 


0.31 


1 1/2 


0.24 


0.45 



Table 2: Woods-Saxon and spin-orbit potential parameters for the continuum final states. 



V 


ro 


a 


v so 




(MeV) 


(fm) 


(fm) 


(MeV) 


(fm) 


-39.8 


1.27 


0.75 


7.07 


0.75 



To describe the valence neutron in the 10 Li continuum we assume that 
the single neutron hamiltonian with respect to 9 Li has the form 

h = t + U + iW (10) 

where t is the kinetic energy and 

U(r)=V WS + SV (11) 

is the real part of the neutron-core interaction. In this paper the imaginary 
part is taken equal to zero. Vws is again a Woods-Saxon potential plus 
spin-orbit whose parameters are given in Table [2} and SV is a correction [17] : 

SV(r) = 16ae 2(r - R)/ 7(l + e {r ~ R)la Y (12) 
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which originates from particle-vibration couplings. They are important for 
low energy states but can be neglected at higher energies. The above form 
is suggested by a calculation of such couplings using Bohr and Mottelson 
collective model of the transition amplitudes between zero and one phonon 
states. Therefore our structure model is not a simple single-particle in a 
potential model but contains in it the full complexity of single-particle vs. 
collective couplings. A more realistic treatment would require the description 
of both bound and unbound states in a three-body model such as in Refs. 
[36] and [H]. 

The continuum energies can be adjusted by varying the parameter a in 
the potential. By changing the strength a of the SV potential in Eq. (1121) we 
could make also the 2s-state just bound near threshold and see what would 
be the effect on the continuum spectrum (cf.[U[8]). As possible final states we 
have considered only the s, p and d partial waves calculated in the potential 
of Table [2] plus Eq. (fT2l , according to Ref . [17] , with different values of 
the strength a. Table [3] gives the energies and widths of the lpi/2 and ld 5 / 2 
states and the corresponding values of a. The widths are obtained from the 
phase shift variation near resonance energy, according to d6j/def\ Erea = 2/Tj, 
once that the resonance energy is fixed [55]. Notice that the values of a for 
the s and d states are very similar, therefore these two states are basically 
obtained in the same potential. 

Table 3: Scattering length of the 2s continuum state, energies and widths of the p- and 
d-resonances in 10 Li and corresponding strength parameters for the SY potential. 





c res 

McV 


r, 

MeV 


a s 
fm -1 


a 
MeV 


2S1/2 
lPl/2 

ld 5 / 2 


0.63 
1.55 


0.35 
0.18 


-17.2 


-10.0 
3.3 
-9.8 



The delta-interaction strength used in Eq. ([6]), is t>2=-8625 MeV fm 3 . It 
has been obtained by imposing that this interaction gives the same volume 
integral as a n- 12 C Woods-Saxon potential of strength -50.5MeV, radius 2.9 
fm and diffuseness 0.75 fm. 
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• Exp. data from [2] 




E f (MeV) 



Figure 2: n- 9 Li relative energy spectrum, for the reaction n Li+ 12 C — > n+ 9 Li+X at 
264 A. MeV. Only the contributions from an s and p initial state with experimental spec- 
troscopic factors [4] C 2 S= 0.3f and 0.45 respectively are included. The thin solid curve 
is the total calculated result. The thick solid curve curve is after convolution with the 
experimental resolution function. The thin dashed curve is the calculation without the 
d-resonance while the thick dashed curve is the same calculation after convolution. The 
symbols with error bars are the experimental points from [2] . Calculations are normalised 
to the data. 

3.2 Results 

Results obtained with the model outlined in Sees. 2 and 3.1 will now be 
discussed. We consider the knockout of a single neutron from a bound state 
in a potential, similarly to the previous calculation for n Be, 14 Be and 14 B pQ. 
The reaction corresponds to a neutron initially bound in n Li which is then 
excited into an unbound state of 10 Li, assuming that another nucleon has 
been emitted and stripped by the target, thus not detected in coincidence 
with the core. One of the results of Ref. p], as far as the reaction model is 
concerned, was to show that projectile fragmentation invariant mass spectra 
depend very weakly on the incident energy and on the neutron initial binding 
energy. Related to this was the investigation of the validity of the sudden 
approximation and the accuracy necessary in calculating the phase shifts. 

In Fig. 2 the symbols with error bars are the experimental points from [2]. 
The thin full curve is the calculated total spectrum, sum of the contributions 
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e f (MeV) 



Fi gure 3: The individual transitions, bound to unbound from the s and p initial state 
components, taken with unit spectroscopic factors, to the s,p and d unbound state as 
indicated. The thick dashed line is the sum of all transitions, including the spectroscopic 
factors for the initial state components while the thick solid line is does not include tran- 
sitions to the d-rcsonancc. 



from the s and p components of the initial wave function, with the parameters 
of Table [T] and the experimental spectroscopic factors from Ref. [4]. Using 
the theoretical spectroscopic factors from Refs. [16], p2] gives very small 
differences in the shape of the total spectrum which however become unno- 
ticeable after convolution with the experimental resolution function. The s, 
p and d final continuum states have the parameters of Table [31 Notice that 
these parameters are perfectly consistent with those extracted from the data 
and given by Table 4 of Ref. [2] . The thick solid curve is the same calcula- 
tion after convolution with the experimental resolution function. The dashed 
curve is the calculation without the d-resonance after convolution while the 
thin dashed curve is the bare calculation without the d-resonance contribu- 
tion. The agreement between the thick solid curve and the data is quite 
good and it shows the importance of including the d-resonance to reproduce 
the tail of the experimental spectrum. The Coulomb breakup spectrum from 
the s component of the initial state was also calculated but since it gives a 
contribution of less than 5% to the peak of the spectrum, it has been omitted 
from the figures. 
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Time dependent a=-17.2fm 

— ■ Eq.(20) of Ref.[3] xO.094 



I , 1 , 1 , 1 , 1 , 1 

0.5 1 1.5 2 2.5 3 

e f (MeV) 

Figure 4: Comparison for the s to s transition of our calculation (solid curve) with that 
according to the sudden formula Eq. (20) of [3|. Both calculations use the exact phase 
shifts. 

The individual transitions, bound to unbound from each initial state com- 
ponent to each possible unbound state as indicated, are shown in Fig. 3. We 
give also in Fig. 4 the comparison for the s to s transition of our calculation 
with the calculation according to the sudden formula Eq. (20) of [3j which 
we discussed at length and compared to our model in Ref. pQ. Both calcula- 
tions were done with the exact phase shift, obtained by solving numerically 
the Schrodinger equation with the potential given by Eqs. ( fTTlfT2l . Contrary 
to what it was found in Ref. P for 13 Be, we find here that in the case of 
the 10 Li virtual state, the scattering length is large enough, i.e. the state is 
close enough to threshold, to justify the use of the sudden approximation. 
Furthermore we have checked that the effective range formula (cf. discussion 
after Eq. (41) of Ref. [1]) gives in this case a very good fit to the exact phase 
shift. In fact, using it in Eq. (20) of [3] gives a curve almost indistinguishable 
from the dashed line in Fig. 4. The parameters obtained from the fit are: 
a s =-14.8 fm, r e = 7.5 fm. Notice that the scattering length value obtained 
instead as a s =-lim^o ^ and given in Table [3] was -17.2 fm. This is the kind 
of sensitivity one can find using the two different prescriptions for a s . On the 
other hand we have checked that the value of r e from the fit, is consistent 
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with the behaviour of our continuum state potential, Eqs. (jll|12j) . Such 
potential becomes indeed negligible, and consistent with zero, for r>r e = 7.5 
fm, as prescribed for the applicability of the effective range theory. 

As already found in Ref. pQ for 11 Be and 13 Be, we can confirm here 
with the 10 Li example that, a part for the s to s transition, the excitation 
of resonances with I > in the continuum of unbound nuclei in projectile 
fragmentation reactions, is a final state effect due to n-core interaction, rather 
than a process in which a bound component of the initial wave function 
becomes suddenly unbound. The results in Fig. 3 show clearly that the 
population of continuum resonance is dominated by the contributions from 
the s-initial state, while the transitions p to p or p to d are not large enough 
to explain the experimental spectrum. This is particularly clear for the d- 
resonance whose presence is necessary to explain the experimental spectrum 
tail but whose strength could not be justified by a d to d transition which 
would have a very small amplitude. Therefore the strength of the continuum 
resonances of a daughter nucleus does not reflect directly the occupation of 
the bound states of same angular momentum in the mother nucleus. This 
is different from the common wisdom on the breakup of two-neutron halo 
nuclei, as discussed for example in the recent Ref. [TJ. 

In the case of 11 Li the two neutrons are in the same state for each compo- 
nent of the initial wave function. Therefore since the p and d wave functions 
have less pronounced tails, the stripping probability of one of the two neu- 
trons, as discussed in Sec. 3 of Ref. [TJ, will naturally diminish the absolute 
value of the peaks due to transitions from these states, with respect to peaks 
due to transition from the bound s-component. This effect is not taken into 
account at the moment in our numerical implementations of the model. On 
the other hand our absolute cross sections should be multiplied by a factor 
two to take into account the fact that the experimental data do not distin- 
guish the two neutrons in the continuum. The absolute value of our cross 
section can be read from the dashed curve in Fig. 3. Taking into account the 
factor two just mentioned, we see that in order to compare to the experimen- 
tal data in Fig. 2, which are given on their absolute scale, we still have to 
renormalise our calculations by a factor two. Considering the incertitude in 
the value of the neutron-target delta-interaction potential and in the strong 
absorption radius, we can consider our estimates quite reasonable. 
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4 Conclusions and outlook 



In this paper we have studied the energy spectrum of 10 Li as obtained from 
the fragmentation of n Li by applying a model [1] for one neutron excitations 
from a bound initial state to an unbound resonant state in the neutron-core 
low energy continuum. The model, based on a time dependent perturbation 
theory amplitude, was previously used pQ to study 13 Be and the breakup of 
n Be and it proved to be reliable. The same conclusion can be drawn here 
after comparing our calculations with the new data from [2] . 

The initial state spectroscopic factors in n Li are quite well known exper- 
imentally, therefore the absolute values of our cross sections have also been 
checked, besides the shape of the n-core relative energy spectrum. Due to 
the closeness to threshold of the s virtual state we have seen that the sudden 
formula used in Ref. [3] and the effective range approximation to the phase 
shift, are both very well justified for 10 Li, contrary to what it was previously 
PP found for 13 Be. Finally we have found that, in agreement with the in- 
terpretation given by the authors of [2J, their recent data provide a clear 
evidence on the excitation of a d resonance around 1.5 MeV. Such a reso- 
nance does not play much role in the composition of the n Li ground state 
but it is an important building bloc of its excited states. 
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